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We study two distinct multi-level atomic models in which one transition is coupled to a Markovian 
reservoir, while another linked transition is coupled to a non-Markovian reservoir. We show that by 
choosing appropriately the density of modes of the non-Markovian reservoir the spontaneous emission 
to the Markovian reservoir is greatly altered. The existence of 'dark lines' in the spontaneous 
emission spectrum in the Markovian reservoir due to the coupling to specific density of modes of 
the non-Markovian reservoir is also predicted. 
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Q\ . I. INTRODUCTION 

T-H ■ 

f^ i' It is now well understood that spontaneous emission of a quantum system depends crucially on the nature of the 

1^ \ reservoir with which the system interacts. Spontaneous emission can be modified in a "tailored" manner by changing 

^^ • the density of modes of the reservoir m . An effective method to achieve this is to place atoms in waveguides [0-0 , 

t~^ _ microcavities [pH10| or photonic band gap materials |ll|-Q, where the density of modes differs substantially from 

CN ■ that of the free space vacuum. The latter has also attracted attention for its potential for modifying the absorption 

and dispersion properties of a system Q . 

\^ , In the above studies the typical scheme involves the interaction of a two-level atom with a reservoir with modified 

J^ ' density of modes. Spontaneous emission of this two- level system differs considerably from the free space result and the 

r^ \ usual Weisskopf-Wigner exponential decay [^ is violated. In specific cases complete inhibition of spontaneous decay 

f^ ■ has even been predicted |^,|3[. In addition to studies of two-level atoms, there are also investigations of multi- level 

\^ \ atoms, where the spontaneous emission from an atomic transition in the modified reservoir infiuences the spontaneous 

^^ . emission of another atomic transition that interacts with the normal free space vacuum |T2,nq,0] . It is schemes such 

CT^ '■ as these that concerns us in this article. Specifically, we study two distinct multi-level atomic schemes where one 

5^ : transition is considered to decay spontaneously in a modified reservoir and the other decays to a normal free space 

r^ ' vacuum. The main interest here is the existence of 'dark lines' (complete quenching of spontaneous emission for 

Oh, specific vacuum modes) in the spontaneous emission spectrum of the free space transition due to the coupling of the 

^.j ' other transition with the modified reservoir. We note that dark lines in spontaneous emission have been predicted in 

Ch \ several laser driven schemes J21-24 and have been experimentally observed |E^. In this article we employ specific 



models for the modified reservoir density of modes, such as, for example, that which results from an isotropic photonic 
band gap model with (or without) defects. We show that dark lines appear in the spontaneous emission spectrum as 
a consequence of the structure of the modified reservoir. 

This article is organized as follows: in the next section we consider a three-level, A-type atomic system, with one 

k^ , transition spontaneously decaying in the normal free space vacuum and the other transition decaying in a modified 

5_] ' reservoir. We study the spontaneous emission spectrum of the free space transition for four different density of modes 

Cu ' of the modified reservoir and show that dark lines can occur in the spectrum due to the structure of the modified 

reservoir. We note that John and Quang ||12| have studied the same atomic system as us, using specifically the 

appropriate density of modes obtained near the edge of an isotropic photonic band gap [see Eq. (^]. However, in 

their study they focussed on the phenomenon of dynamical splitting in the spectrum and not on the existence dark 

lines, which is the main phenomenon discussed in our article. Zeros and splittings in spectra should not be confused, 

although they relate to each other in the limit of strong coupling. In section III we consider a laser-driven extension 

of the previous system and show that dark lines can occur in this system, too. In this case the dark lines originate 

from either laser-induced or modified vacuum-induced mechanisms. Finally, we summarize our findings in section IV. 

II. FIRST CASE: A-TYPE SCHEME 

We begin with the study of the A-type scheme, shown in Fig. 0(a). This system is similar to that used by Lewenstcin 



et al [H| and by John and Quang |12|. The atom is assumed to be initially in state |2). The transition |2) ^^ |1) is 
taken to be near resonant with a modified reservoir (this will be later referred to as the non-Markovian reservoir), 
while the transition |2) <-^ |0) is assumed to be occurring in free space (this will be later referred to as the Markovian 
reservoir). The spectrum of this latter transition is of central interest in this article. The Hamiltonian which describes 



the dynamics of this system, in the interaction picture and the rotating wave approximation (RWA) , is given by (we 
use units such that ?i = 1), 



^ = E 



9\e 



— i{uix—'^2a)t 



|2)(0|aA+5].9.e-^("~-'^-)*|2)(l|a, + H.c. 



(1) 



Here, g^ denotes the coupUng of the atom with the modified vacuum modes (k) and g\ denotes the couphng of the 
atom with the free space vacuum modes (A). Both couphng strengths are taken to be real. The energy separations of 



the states are denoted by coi 



LUj and lUk i'-^x) is the energy of the k (A)-th reservoir mode. 



The description of the system is given using a probabiUty ampUtude approach. We proceed by expanding the wave 
function of the system, at a specific time i, in terms of the 'bare' state vectors such that 



mt)) = him, {0}) + E ^^ wiO' {^}) + E ^^ wii' {'*» ■ 

A K. 

Substituting Eqs. (|l|) and (||) into the time-dependent Schrodinger equation we obtain 

X K, 

ib,{t) = g^b2{t)e'^'^^-^^'^' . 



(2) 



(3) 

(4) 

(5) 



We proceed by performing a formal time integration of Eqs. (Q) and (g) and substitute the result into Eq. (^) to 
obtain the integro-differential equation 



Ut) 



dt'b2{t')^gle 



2 -i(wA-'^2o)(t-t') 



dt'b2{t')Y,9lt 
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Because the reservoir with modes A is assumed to be Markovian, we can apply the usual Weisskopf-Wigner result |2C 
and obtain 



■,2p-J('^A-W20)(t-t') _ l-^d — 



E 5^6-'^"'^ -"^"^^*-*'^ = ^5{t - t') . 



(7) 



Note that the principal value term associated with the Lamb shift which should accompany the decay rate has been 
omitted in Eq. (|^). This does not affect our the results, as we can assume that the Lamb shift is incorporated into the 
definition of our state energies. For the second summation in Eq. (o), the one associated with the modified reservoir 
modes, the above result is not applicable as the density of modes of this reservoir is assumed to vary much quicker 
than that of free space. To tackle this problem, we define the following kernel 



K{t - *') = E .9^e-'(— -)(*-*') « 5^/2 j dup{oj)e 



i{uj—uj2i){t — t') 



(8) 



which is calculated using the appropriate density of modes p{i^) of the modified reservoir. In Eq. (0), g denotes the 
coupling constant of the atom to the non-Markovian reservoir. Using Eqs. (0) and (g) into Eq. (|^) we obtain 



Ht) 



iUt) 



dt'b2{t')K{t~t') 



(9) 



The long time spontaneous emission spectrum in the Markovian reservoir is given by S{5x) oc \b\{t —>■ oo)p, with 
S\ = ijj\ — ti;20 l|2l|,^|. We calculate 6a (t —^ oo) with the use of the Laplace transform Ig^l of the equations of 
motion. Using Eq. (^) and the final value theorem P6[ we obtain the spontaneous emission spectrum as S{6x) oc 
7I lims-»-i(5A 52(s)P, where B2{s) is the Laplace transform of the atomic amplitudes b2{t) and s is the Laplace variable. 
This in turn, with the help of Eq. (pf) , reduces to 



S{5x) (X 



7 



-iSx + 7/2 + k{s 



-i5\ 



(10) 



where K{s) is the Laplace transform of K{t), which yields 

K{s) = g3/2 f ^^ 



p{uj) 



s + i{uj — UJ21) 



(11) 



Therefore, in order to calculate the spontaneous emission spectrum in the Markovian reservoir we need to calculate 
K{s). This will be done for different models of the density of modes p{uj) of the non-Markovian reservoir. 

We begin by considering the non-Markovian reservoir to be that obtained near the edge of an isotropic photonic 
band gap model |l|,|ll. Then, 



p{uj) 



1 



-.e{uj-ujg), 



TT ^LO - LOg 

where ujg is the gap frequency and 8 is the Heaviside step function. In this case, Eq. (pi]) leads to 



K{s) = 



,3/2 



V** - ^9 ' 



with 5g = LOg — L021- The spontaneous emission spectrum then reads 



S{5x) oc 7 





y/5\ - 5g 




(- 
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5g + g^/^ 



(12) 



(13) 



(14) 



Obviously, the spectrum exhibits a zero (i.e. predicts the existence of a dark line), if 5\ = Sg. This is purely an 
effect of the above density of states, and the non-Markovian character of the reservoir. In the case of a Markovian 
reservoir the spontaneous emission spectrum would obtain the well-known Lorentzian profile and no dark line would 
appear in the spectrum. The behaviour of the spectrum is shown in Fig. for different values of the detuning from 
the threshold. The spectrum has two well-separated peaks and the dark line appears at the predicted value. We note 
that a spectrum of this form has also been derived by John and Quang |l3|; however, in their study they did not 
mention the existence of the dark line (i.e. the zero in the spectrum and its spectral position), which is of central 
importance in this article. 

The density of modes of Eq. (03) is a special case of a more general family of density of modes those given by 



p{uj) = ^ — 



-UJo 



-Q{uj~ujg), 



(15) 



where e is usually referred to as the smoothing parameter. Such a density of modes has been used in studies of atoms 
in waveguides JEIP] microcavities P,p| and photonic band gap materials M] . Eq. (O) is recovered by taking the limit 
e ^ in Eq. (|15[). For the above density of modes, Eq. (p^) we obtain 



K{s) 
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(16) 
(17) 



So, in this "smoothed" case, no dark line appears in the spectrum. This is shown in Fig. H, where the zero disappears 
from the spectrum. The only case for which spontaneous emission spectrum will give zero is that of e = and S\ — Sg 
which, of course, reduces to the previous result of Eq. (14). 

The above two density of modes describe "pure" materials (i.e. materials with no defects). In reality, defects exist 
in waveguides, microcavities or photonic band gap materials which exhibit gaps in their density of modes. These 
defects lead to narrow-linewidth, well-localized modes in the gaps of the above structures. Their density of modes can 
be described by either a delta function pl| ] or a narrow Lorentzian 1 13[ . We will now combine the density of modes 
given by Eq. (p^ ) and that given by a defect structure. In the case that we choose a delta function density of modes 
for the defect structure we obtain 
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(18) 
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where 5c = oJq — ^^-"21 is the defect mode-atom detuning (defect mode at frequency uoc) and gi is the defect mode-atom 
couphng constant. In this case the spectrum exhibits two dark hnes one at 5\ = 5g and another at 5\ = 5c- In Fig. 
^ we present the spontaneous emission spectrum described by Eq. (jl^) for the same parameters that used in Fig. ||. 
The spectrum obtains now a very pronounced third peak and two dark hnes at the predicted values. For a Lorentzian 
density of modes of the defect structure we get 
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(20) 
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with 7c being the width of the Lorentzian describing the defect mode. As in the case of the density of modes given 
by Eq. ([12|), in this case too, there is a single dark line in the spectrum at Sx = Sg. However, the behaviour of the 
spectrum is quite different in this case compared to that shown in Fig. ^ as can be seen in Fig. ||. 



III. SECOND CASE: LASER-DRIVEN SCHEME 



In this section we turn to the study of the laser-driven system shown in Fig. |l|(b). This system is composed of a 
ground state |3) which is coupled by a laser field to the excited state |2). State |2) can spontaneously couple to either 
state |0) via interaction with a Markovian reservoir, or to state |1) via interaction with a non-Markovian reservoir, as 
in the previous section. The spontaneous emission spectrum in the Markovian reservoir is our main concern in this 
section, too. The Hamiltonian of this system, written in the interaction picture and under the RWA reads 



H = n\3){2\e'^'+J2 



gxe 



— i{uJx~^2o)t 



\2){Q\ax 



jKe 



— i{u)K—u)-2l)t 



|2)(l|a« + H.c.. 



(22) 



Here, fi is the Rabi frequency (assumed real) and 5 — u — 0^23 the laser detuning, with u being the laser field angular 
frequency. For the description of this system too we use the probability amplitude approach and expand the wave 
function of the system as 



|V'(i))=63(0e^^'|3,{0})+62(t)|2,{0}) 



E 

A 



6A(t)|0,{A})+^6«(t)|l,{4) 



(23) 



We substitute Eqs. ( p2D and (|23| ) into the time-dependent Schrodinger equation and apply the elimination of the 
vacuum modes outlined in the previous section to obtain 



ih{t)^5h^{t) + nb2{t), 



ib2{t) = nb^it) - ii-b2{t) - ^ / dt'b2{t')K{t ~ t') 
2 Jo 
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(25) 

(26) 
(27) 



The long time spontaneous emission spectrum in the Markovian reservoir is given by S{5x) c>c \bx{t — > oo)p, and is 
calculated in a closed form with the use of the Laplace transform and the final value theorem Igg as 



S{5x) ex 7 



(5a- (5)62(0) + f^fe3(0) 



{6x-5)[6x+ill2 + iK(i 



-iSx)] - r!2 
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The spectrum depends, in this case too, from the Lap lace transform of the kernel, which in turn depends on the 
density of modes of the modified reservoir [see Eq. (|l^)] . We use the four models of the density of modes described 
in the previous section to calculate the above spectrum. 



In the case of an isotropic photonic band gap material, with the use of Eq. (13) we find, 
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Let us assume first that the atom starts from its ground state, i.e. |&3(0)p = 1, |62(0)p = 0. Then, the spectrum 
obtains a dark hne at 6\ = (5g, except in the case that 6 = 6g, so in the case that the laser detuning becomes equal to 
the detuning from the band edge no dark line occurs, if the system is initially in state |3). We note that in the case 
when the transition 1 2) — > 1 1) occurs in a Markovian reservoir, then no dark line appears in the spectrum |2l| , |24| . If now 
the atom starts in a superposition of the ground (|3)) and excited states (|2)) with real expansion cocfhcients [foj(O)] 
then two dark lines can appear in the spectrum. The first dark line appears at S\ = (5 — i763(0)/fe2(0) and is attributed 
to the laser-atom interaction |24|. The second dark line occurs at S\ — 5g and is attributed to the interaction with 



the non-Markovian reservoir. In this case too, the second dark line disappears ii 5 = 5g. The behaviour of the above 
spectrum for the case that the system is initially in state |2) is shown in Fig. H. We have chosen 5g ^ 5 therefore, 
as predicted, two dark lines appears in the spectrum. We note the similarity of this spectrum and that of Fig. IJ. 
This should be expected as the delta function density of modes (used in Fig. 0) represents a pure Jaynes-Cummings 
interaction pTI . 

If the more general, smoothed density of states of Eq. ( p^ ) is used (with e 7^ 0), the spontaneous emission spectrum 
is given by 
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Then, in the case that the atom starts in the ground state no dark line exists in the spectrum. However, if the atom 
starts in a superposition of the ground and excited states (with real expansion coefficients) the laser-induced dark line 
appears at the same frequency as before, i.e. at S\ = S — J763(0)/62(0). This is verified in Fig. ^ where only a single 
dark line appears in the spectrum at dx — S, as 63(0) = 0. 

In the case of an isotropic photonic band gap with a defect with delta function density of modes then the spectrum 
reads 
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and one obtains, in general, two dark lines if the atom starts from the ground state, one at Sx = Sg and the other at 
dx = 5c- Any of these dark lines can disappear if either 5 = 5g, ov 5 = 5c so the spectrum exhibits only one dark line 
in this case. If now the atom starts in a superposition of the ground (|3)) and excited states (|2)) with real expansion 
coefficients then a third dark line can appear in the spectrum aX 5x — 5 ~ r2&3(0) 762(0). The spectrum for the case 
that the atom is initially in the excited state |2) (as in the previous cases), is displayed in Fig. B. A very rich behaviour 
of the spectrum is find. The spectrum now has three dark lines at the predicted values and four different peaks can 
be seen. 

Finally, in the case of an isotropic photonic band gap with a defect with Lorentzian density of modes the spectrum 
gets the form, 
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The dark lines in the spectrum appear at the same frequencies as those of the simple isotropic photonic band gap 
model, discussed above. However, the shape of the spectrum in this case differs from that of Fig. ra as it is shown in 
Fig. |. 



IV. SUMMARY 



In this article we have investigated the spontaneous emission properties of two distinct atomic models with one 
transition coupling to a Markovian reservoir while another transition coupling to a non-Markovian reservoir. Of 
specific interest to us were the existence of dark lines in the Markovian spontaneous emission spectrum. We have 
shown that dark lines can occur if the non-Markovian reservoir is described by certain densities of modes. In the case 
of the laser-driven scheme of Fig. 0(b) laser-induced dark lines can co-exist with non-Markovian reservoir-induced dark 
lines. Overall, spontaneous emission in the Markovian transition can be efficiently controlled (and even suppressed) 
by appropriately engineering the density of modes of the non-Markovian reservoir. 
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FIG. 1. The two systems under consideration. In (a) the 
thick dashed line denotes the coupling to the modified reser- 
voir and the thin dashed line denotes the coupling to the 
Markovian reservoir. The same hold in (b), and in addition 
the solid line denotes the coupling by the laser field. 




Detuning 
FIG. 2. The spontaneous emission spectrum S{5x) (in ar- 
bitrary units) given by Eq. ( |l4| ) for parameters g — 1, and 
5g = (dotted curve); 5g = 1 (dashed curve); Sg — —1 (full 
curve). All parameters are in units of 7. 




Detuning 
FIG. 3. The spontaneous emission spectrum S{S\) given 
by Eq. ( [l7| ) for parameters g — 1, e = 0.3, and 5g — (dotted 



s); S, 



g — 1 (dashed curve 



y,Sg 



-1 (full curve) 




Detuning 
FIG. 4. The spontaneous emission spectrum S{S\) given 
by Eq. (hol) for parameters g — 1, gi — 1, 5c — —2, and 
Sg = (dotted curve); 5g = 1 (dashed curve); Sg = —1 (fuU 
curve) . 
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Detuning 
FIG. 5. The spontaneous emission spectrum S{5x) given 
by Eq. ( |2l| ) for parameters g = 1, (71 = 1, 7c = 1, <5c = —2, 
and (5g = (dotted curve); 5g = 1 (dashed curve); Sg = —1 
(full curve). 




Detuning 
FIG. 6. The spontaneous emission spectrum S{S\) given 
by Eq. Q for parameters g = 1, D, = 1, S = -1.5, 62(0) = 1, 
and Sg — (dotted curve); 5g — 1 (dashed curve); 5g — —1 
(full curve). All parameters are in units of 7. 
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Detuning 
FIG. 7. The spontaneous emission spectrum S{S\) given 
by Eq. ( po|) for parameters g — 1, e = 0.3, Q, = 1, S = —1.5, 
^2(0) — 1, and Sg = (dotted curve); Sg — 1 (dashed curve); 
Sg — ~1 (full curve). 
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Detuning 
FIG. 8. The spontaneous emission spectrum S{S\) given 
by Eq. (Bll) for parameters g = 1, gi = 1, Sc = —2.5, fi = 1, 
5 — —1.5, 62(0) = 1, and 5g — (dotted curve); 5g = 1 
(dashed curve); 5g — —1 (full curve). 
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Detuning 
FIG. 9. The spontaneous emission spectrum S{S\) given 
by Eq. ( |32[ ) for parameters g — 1, gi — 1, 'jc = 1, 5c = —2.5, 
Q = 1, 5 = —1.5, 62(0) = 1, and Sg = (dotted curve); Sg — 1 
(dashed curve); Sg = ~1 (full curve). 
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